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Abstract 

I demonstrate that the chiral properties of Domain Wall Fermions (DWF) in the 
large to intermediate lattice spacing regime of QCD, 1 to 2 GeV, are significantly im- 
proved by adding to the action two standard Wilson fermions with supercritical mass 
equal to the negative DWF five dimensional mass. Using quenched DWF simulations 
I show that the eigenvalue spectrum of the transfer matrix Hamiltonian develops a 
substantial gap and that the residual mass decreases appreciatively. Furthermore, I 
confirm that topology changing remains active and that the hadron spectrum of the 
added Wilson fermions is above the lattice cutoff and therefore is irrelevant. I argue 
that this result should also hold for dynamical DWF and furthermore that it should 
improve the chiral properties of related fermion methods. 



1 Introduction 



Domain Wall Fermions (DWF) and other closely related methods provide the most faithful 
lattice regularization of QCD with unprecedented chiral symmetry and topological proper- 
ties. These methods have produced an impressive wealth of results and have brought about 
a new and brave era for lattice gauge theory. Several proposals for improvement have also 
been demonstrated. The reader is referred to the reviews [TJ|21IIIlSII3ini[ZllBllH], related 
works P0]> and references therein for more details. Here I will focus on DWF and QCD. 
The method I will described should apply to all closely related lattice fermions as well as to 
other gauge theories similar to QCD (such as = 1 Super Yang-Mills |22| ) . The method 
was first proposed in [I]. 

Using the new generation supercomputers and numerical simulation algorithms one can 
now simulate dynamical QCD at zero and finite temperature in the strong to intermediate 
coupling regime with lattice spacing a in the region 1 GeV < a -1 < 2 GeV and for 
number of lattice points L s along the fifth dimension in the neighborhood of L s rj 20 to 100 
depending on the lattice spacing a. These L s values are now certainly within the reach 
of the latest supercomputers. However, one obviously would like to do better. A method 
that would improve DWF so that smaller values of L s are needed would be very welcomed. 
And not only for the obvious benefits in computing time but also as a matter of theoretical 
principle. The reasons behind the need for rather large values of L s are as involved as 
they are interesting. But they are also a lattice artifact. It should be possible to remove it 
without spoiling the important features of DWF. 

In this work I demonstrate that the chiral properties of DWF in the large to intermediate 
lattice spacing regime 1 GeV < a -1 < 2 GeV are significantly improved by adding to the 
action two standard Wilson fermions with mass equal to the negative DWF five dimensional 
mass mo (Wilson supercritical region). These fermions produce a fermion determinant that 
is identically zero when the fifth dimension transfer matrix Hamiltonian if 4 (mo) has a zero. 
Therefore, small eigenvalues of if 4 (m ) ought to be suppressed in the simulation. And it 
is precisely these eigenvalues that are a lattice artifact and are directly responsible for 
the large L s requirement. Suppressing them allows one to achieve same chiral symmetry 
restoration for smaller values of L s . Indeed, I show that the eigenvalue spectrum of #4(7710) 
develops a substantial gap and that the residual mass m res , quantifying the remnants of 
chiral symmetry breaking, decreases appreciatively. I use this key property to name this 
fermion regulator as Gap Domain Wall Fermions (GDWF). For other DWF related fermions 
(for example overlap-Neuberger fermions) I argue that this method of adding gap fermions 
(GF) should also improve their chiral properties. 

In order to make sure there are no adverse "side-effects", I confirm that local topol- 
ogy changing (measured using the overlap index method) remains active; instantons and 
anti-instantons with size larger than a lattice spacing are active in the gauge field configu- 
rations. Also, I confirm that the hadron spectrum of the added Wilson fermions is above 
the lattice cutoff and therefore is irrelevant. Furthermore, I check that parity is not broken 
and therefore the added Wilson fermions are well outside the Aoki phase. All numerical 
simulations in this work were performed on 512- node IBM BlueGene/L supercomputers. 

This paper is organized as follows: In section |21 DWF and the difficulties with chiral 
restoration are briefly reviewed. In section El the Gap Domain Wall Fermion (GDWF) 
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method is presented. In section numerical results using GDWF in the quenched approx- 
imation on 16 3 x 32 lattices are given for three lattice spacings a ~ 1.0, 1.4 and 2.0 GeV , 
and are compared with the quenched approximation of standard DWF at the same lattice 
spacings. Also, in this section I argue that these results will hold for dynamical GDWF as 
well. In section an example of a quenched two-flavor QCD simulation with GDWF at 
the rather large lattice spacing a -1 pa 1.4 GeV is presented. It is shown that the pion mass 
acquires physical values already for an L s as small as 24 with m res contributing about 10% 
to the bare quark mass. In section El I discuss issues that relate to topology. In section 
the ease of implementing GDWF, as well as the small GDWF computational cost added to 
the cost of a DWF simulation is discussed. In section |S1 some curious theoretical thoughts 
brought up by GDWF are pondered upon. In section El I discuss the applicability of the 
method to overlap and Neuberger fermions and furthermore to the Af = 1 Super Yang-Mills 
theory. Also, I present open issues for further exploration. The conclusions are given in 
section [TUJ 

2 DWF and the chiral restoration problem 

Lattice DWF ^TJ ^3 El EH] are defined in five dimensions. The fifth dimension has 
L s lattice sites and the five-dimensional fermion has positive bare mass m (domain wall 
height). The five dimensional Dirac operator Dp employs free boundary conditions at the 
edges of the fifth dimension (walls). As a result the plus chirality fermionic components are 
localized on one wall while the minus chirality components are localized on the other. The 
two chiralities are explicitly mixed with a mass parameter m f . The gauge fields are defined 
in four dimensions only. They are the same along the fifth dimension and have no fifth 
component. This allows for a definition of a transfer matrix T along the fifth direction that 
is the same in all "slices" along that direction. The product of the transfer matrices along 
the fifth direction is therefore T Ls . The single particle Hamiltonian if 4 (mo) associated 
with this transfer matrix is then also independent of the fifth dimension. It is defined in 
four dimensions and, for the case where the fifth dimension is continuous, one can show 
that if 4 (m ) = 75-^ (— m ) where p w (— m ) is the standard Wilson fermion Dirac operator 
with mass —mo- When the fifth dimension is not continuous the Hamiltonian has a more 
complicated form, but one can show that it has the same zero eigenvalues as H^mo). 
The five dimensional Dirac operator is given by the equation: 



Dp(x, s; x' , s'; mo, m/) = 5(s — s')lp{x, x'\ m ) + Jp^{s, s'\ mf)8{x — x'\ mf) 



(1) 



p(x,x';m ) 



J2{(l + l»)U,(x)6(x + f,-x') + (l 



lv)Ul{x')5{x' + ft -x) 



(4 — m )5(x — x') 



(2) 



p-(s, s';m f ) 



P R 5(1 - s') - m f P L 5(L s - 1 - s') - 5(0 - s') s = 

P R 6(s + l-s') + P L 8{s -l-s')-8(s-s') < s < L s - 1 

-m f P R 5{0 - s') + P L 5(L S - 2 - s') - 5{L S - 1 - s') s = L s - 1 



(3) 
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where mo is the 5-dimensional mass representing the "height" of the Domain Wall and m/ 
is the explicit bare quark mass. The gamma matrices 7 are taken in the chiral basis. 
The standard Wilson fermion Dirac operator is: 



The localization of the two chiral components on the opposite walls is exponentially 
good. The slowest decay coefficient is proportional to the smallest, in absolute value, nega- 
tive eigenvalue of #4 (mo). For infinite L s (overlap fermions) the two chiralities completely 
decouple provided that i7 4 (m ) does not have eigenvalues that are exactly zero. That pos- 
sibility is of measure zero and is therefore of no concern. Nevertheless, at finite L s , where 
simulations are performed, the two chiralities will mix and break chiral symmetry. This 
mixing is of a similar nature as the one produced by a mass term. It is possible to calculate 
this "effective" mass (usually called residual mass, m res ) and use it to quantify the quality 
of the DWF regulator. Clearly at finite L s one would like H^tuq) to have a substantial 
gap which in turn would result to a rapidly decreasing increased. 

Obviously the spectrum of H^fi) at various values fj, is of particular importance. A 
small review is given below. It will serve to explain GDWF. The reader is referred to the 
original work [T2] for more details. For any gauge field configuration H^(fi) has the same 
number of positive (n+) and negative (n_) eigenvalues for fi < 0. However, as \i is increased 
above zero some eigenvalue of H±{p) may cross zero and change sign. Then n + — n_ would 
not be zero just after the crossing occurs. It has been shown that the number and direction 
of crossings is directly related to the number of instantons and anti-instantons present in 
the gauge configuration [T21 EH! and that n + — n_ is in fact equal, in a statistical sense, to 
the net (global) topological charge of the gauge field configuration j~fj""]. The Atiyia-Singer 
index theorem is realized on the lattice in a statistical sense. These are rather remarkable 
properties. 

A very nice way to look at the spectrum of H 4 (fj,) is to plot the eigenvalues of H^/i) as a 
function of /1. This is an eigenvalue flow diagram [1211131221 (see Figure 4). Instantons that 
are larger than the lattice spacing are of course of physical interest and it has been shown 
[2*3] that they produce crossings in rather localized neighborhoods of \i that correspond 
to the edges of the standard Brillouin zones. For a single flavor DWF one picks ji = itlq 
in-between the first and second set of crossings. Since this is a finite range no fine-tuning is 
required. In the continuum limit the range for one flavor extends from to 2 and the width 
of the neighborhoods where crossings occur tends to zero. The location and width of the 
crossing neighborhoods is renormalized from their continuum values as the coupling is made 
stronger. For example, for lattice spacing a -1 ~ 1.4 GeV , the first set of crossings occur in 
the neighborhood of /i m j n ~ 0.9 and the second in the neighborhood of /i max ~ 2.2. Their 
width is approximately 0.2. However, small instantons of the size of the lattice spacing are 
generated/destroyed because they can "come up/fall in" through the discrete lattice. This 
generates additional crossings throughout the /i m i n , /i max region. 

In a numerical simulation at small lattice spacing (weak coupling) there are few to no 
such small instanton crossings. The simulation is performed at a fi = mo in the middle of 




(5) 
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the relevant range. At that value, since there are no crossings, the eigenvalue gap is large 
and therefore the localization on the walls is good. As a result, the two chiralities mix 
very weakly and break chiral symmetry minimally (as a result m res decreases rapidly with 
increasing L s ). However, at large lattice spacings (strong coupling) the small instantons 
generate crossings across the whole range and therefore also close to mo- As a result the 
eigenvalue gap becomes very small. To be sure the probability of an eigenvalue crossing 
exactly at m is zero. There is always a gap but it may be very small at large lattice spacings. 
The challenge is to suppress the crossings due to the lattice spacing size instantons, which 
are obviously an artifact of the lattice "discreteness" , without destroying the crossings due 
to the all-important physical instantons with size of many lattice spacings. 

3 Gap Domain Wall Fermions 

A method that would produce a sizeable gap in the spectrum of H^rrio) for lattice spacings 
in the 1 GeV ^ a -1 ^ 2 GeV region would be of great use. Here I describe such a method. 
Gap Domain Wall Fermions are similar to standard Domain Wall Fermions but they have 
the desired gap. 

Since if 4 (m ) = ^^D w { — m ) where D w (—m ) is the standard Wilson fermion Dirac 
matrix one can induce a larger gap by adding to the theory standard dynamical Wilson 
fermions with mass —mo. Here I add two flavors. When integrated out these fermions 
contribute a factor of det 2 [D w (—m )} = det 2 [H^(mo)] to the Boltzman weight. Gauge field 
configurations for which if 4 (mo) has small eigenvalues will be suppressed by this Boltzman 
weight and therefore they will be sampled very infrequently. In particular, any gauge field 
configuration for which if 4 (mo) has a zero eigenvalue is explicitly excluded (not to mention 
that the set of such configurations is of measure zero). A gap at m is ensured. But more to 
the point is the fact that this Boltzman weight "repels" gauge configurations for which the 
gap at mo is small. One can expect a substantially larger gap even for strong couplings. In 
section 01 I present numerical results that demonstrate that indeed this is the case. Notice 
that I have not added any extra parameters since the Wilson fermions have mass equal to 
the negative five dimensional mass which is already a parameter of the theory. 

The Wilson fermions that I added to the theory have mass —mo with mo somewhere 
in the middle of the crossings region [/i m m, /^max]- I have chosen m = 1.9 which is a good 
choice for the whole range of lattice spacings of interest. Such a mass is in the supercritical 
region of Wilson fermion masses and is very heavy. The hadron spectrum, including the 
pions, of these two flavors of Wilson fermions, should be above the cutoff. In that case their 
contribution to the low energy physics of the theory is irrelevant. Indeed, in section HJ I 
present numerical results that demonstrate that their hadron spectrum is above the cutoff. 
Furthermore, I show that parity is not broken and therefore the added Wilson fermions are 
well outside the Aoki phase. 

As mentioned earlier, it is important that crossings due to the all-important physical 
instantons with size of many lattice spacings are present. The added Wilson fermions have 
mass —mo and they suppress the crossings around mo but have little effect further away. 
Because mo is chosen somewhere in the middle of the allowed range, the larger instanton 
crossings are not affected since they occur at the edges of the allowed range. Again, in 
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section I present numerical results that demonstrate that this is the case. 
The GDWF full partition function for QCD is given by: 

Z = J [dU] J [dXdX] J [dVdV] J [d&d<S>]e- s (6) 

Upix) is the gauge field, X(x) is the 4-dimensional Wilson fermion field ^f(x, s) is the 5- 
dimensional fermion field, and $>(x,s) is the 5-dimensional bosonic Pauli Villars (PV) type 
field, x is a coordinate in the 4-dimensional space-time box with extent L along each of the 
directions, \i — 1, 2, 3, 4, and s — 0, 1, . . . , L s — 1, where L s is the size of the fifth direction 
and is taken to be an even number. The action S is given by: 



S = S(/3, L, L s , m ,m f ) = 



S G (U;P,L) 
S W (X,X, U; -mo, L) 
S F (y,y,U,m ,m f ,L,L s ) 
+ S PV {&,$,U;m ,L,L s ) (7) 



+ 
+ 



where: Sq, is the standard plaquette Wilson gauge action with coupling g (/3 = 6/g 2 ) EHl. 
Of course, any pure gauge action can be chosen instead. In particular, one can choose 
a lattice gauge action that most closely resembles the continuum action. For example, a 
popular choice is the Iwasaki action [22]. Because, as mentioned above, DWF maintain a 
connection between their index and the topological charge, it is possible to improve them 
by suppressing small instantons by an appropriate choice of the lattice pure gauge action 
[H El I2H1 EO] • The Iwasaki type actions have been shown to have this property but mostly 
for lattice spacings in the lower end of the interval considered here and below. In any 
case, in order to be able to evaluate GDWF alone I did not want the improvements due 
to these actions to obscure the results and that is why I use the standard Wilson gauge 
action. However, it should be obvious that for production simulations one should use an 
improved gauge action not only because of the possible additional improvement on the 
chiral properties but also because of the closer resemblance to the continuum gauge action. 
Sw is the standard Wilson fermion action [21] with mass — m . Sf + Spy is the standard 
Domain Wall Fermion action with the Pauli Villars regulator, five dimensional mass uiq 
and explicit bare quark mass m f [TH H2l HSl HI] • 



4 Numerical results 

In this section I present numerical results that demonstrate the properties of GDWF. 

Because of limited computational resources I use the "quenched" approximation for the 
DWF sea fermions. However, unlike standard quenched simulations I consider large lattice 
spacings that are of interest to dynamical DWF simulations. I compare results obtained 
from simulations with no dynamical Wilson flavors with results obtained at the same lattice 
spacing with two dynamical Wilson flavors with mass —tuq. I only use one value of mo = 1.9 
throughout this work. In order to compare results I match the lattice spacing between the 
two cases by adjusting (3. Here I achieve a 5% or better matching level at three values of 
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the lattice spacing a -1 ~ 1.0 GeV , 1.4 GeV and 2.0 GeV . Measurements are done using 
the DWF operator at m = 1.9. The space-time volume of all simulations is 16 3 x 32. More 
technical details about the numerical simulations are given in the appendix. 

Therefore, the two cases correspond to standard quenched DWF simulations and to 
quenched GDWF simulations. Again, the lattice spacings I use are large and of interest 
to dynamical simulations. I expect that the effects on chiral symmetry due to the DWF 
quenching are to a large extent taken into account by simulating at these large lattice 
spacings. To be sure the quenched approximation does not allow for the calculation of a 
systematic error. However, if previous results are of any guide I expect that quenching 
may affect mass measurement results at the 10% level. Since the improvements observed 
here are at the one order of magnitude level I do not expect that quenching is affecting the 
results in any significant way. 

The first order of business is to match the lattice spacings between DWF and GDWF. 
I measure the p and nucleon masses at L s = 16 and extrapolate to rrif = 0. The results for 
three values of (3 are shown in Figure 1. These values of (3 were chosen so that the rrif = 
extrapolated p masses match the quenched DWF masses of |27j . 

The matching for the rrif — extrapolated masses is shown in Figure 2. The DWF 
p masses match well (better than 5%). I choose the DWF rrif = extrapolated p mass 
to set the scale. The values of j3 shown correspond to a -1 « 1.0 GeV , 1.4 GeV and 
2.0 GeV . The nucleon masses also match well (better than 5%) which is an indication 
that the method works properly. From now on I will use these values of beta to perform 
comparisons of the chiral properties of DWF and GDWF. The inverse lattice spacing in 
GeV corresponding to the quenched GDWF (3 values is plotted in Figure 3. 

In Figure 4 the ten smallest magnitude eigenvalues of if 4 (mo) are plotted vs. m . 
The eigenvalues are calculated with an accuracy 10 -6 and are measured in m steps of 
0.025. An aggregate of the results from 20 independent configurations (separated by 20 
configurations) is plotted in each plot. The left column is from 0-flavor Wilson simulations 
while the right column is from 2-flavor Wilson with mass — m . Horizontally, the 0-flavor 
and 2-flavor (3 values correspond to the same lattice spacing (Figure 2). From top to bottom 
a -1 « 1.0 GeV , 1.4 GeV and 2.0 GeV . Notice the difference in the y-axis scale for the 
different lattice spacings. The "crosshairs" indicate the mo = 1.9 point. One can clearly 
see that the 2-flavor Wilson fermions generate a substantial gap around mo where none 
existed before even at the large lattice spacing a -1 « 1.0 GeV . 

Furthermore, it is very important to observe in Figure 4 that the 2-flavor Wilson 
fermions generate the gap at a neighborhood of m = 1.9, but allow for a copious amount 
of crossings at the edges of the allowed mo range. As mentioned in section |2l these crossings 
correspond to instantons with size larger than a lattice spacing and are of physical interest. 
Although Figure 4 shows the cumulative results of 20 configurations, by close inspection I 
confirmed that the number of crossings changes from configuration to configuration. This 
indicates instanton, anti-instanton activity. However, with the current resolution I am not 
able to measure the net topological index. This is beyond the scope of this work. For 
further discussion see section |B] 

In order to get a better picture of the small eigenvalue distribution I chose to look at 
the m = 1.9 "cross-section" for the and 2 flavor cases at a" 1 rj 1.4 GeV {(3 = 5.85, 
4.6) using a larger number of eigenvalues and configurations. I measured the 100 smallest 
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magnitude eigenvalues of #4 (mo) at m = 1.9 for 110 independent configurations (separated 
by 20 configurations, i.e. a total of 2200 configurations were generated). The fraction of 
eigenvalues with values between A and \ + d\ is plotted vs. A in Figure 5. Here dX ~ 0.007. 
The difference in the distributions is telling. The 0-flavor Wilson distribution has non-zero 
support at A = and is raising almost linearly with a small slope as A increases. On the 
other hand, the 2-flavor Wilson distribution has zero support between A = and A = 0.015. 
After that it is raising very slowly until A = 0.05. It does not exceed the 0-flavor distribution 
until A gets larger than pa 0.12. Larger statistics are needed to paint a clearer picture. Of 
course it may be that the distribution has zero support only at A = and then it increases 
perhaps quadratically. Nevertheless the point is that for a numerical simulation of QCD 
that typically generates a few thousand configurations the method creates a substantial 
gap. 

As discussed, it is expected that the added 2-flavors of Wilson flavors with supercritical 
mass /i = —1.9 should have a hadron spectrum above the lattice cutoff. This is verified 
in Figure 6. The pion (diamonds), rho (squares) and nucleon (stars) masses vs. a -1 in 
GeV measured using 2-flavor dynamical Wilson fermions (both for the propagator and the 
sea quarks) is shown. The scale is set using the GDWF p mass. The straight line marks 
the cutoff. Masses above that line are above the lattice cutoff. Clearly the Wilson hadron 
spectrum is above the lattice cutoff. Most importantly the pion mass is very heavy and 
above the cutoff. This ensures that the added 2-flavors of Wilson fermions do not affect 
the low energy physics. On the other hand, they clearly affect the physics near the cutoff. 
For example, they renormalize the value of (3. At a" 1 ~ 1.4 GeV , (3 changes from 5.85 to 
4.6. But, of course, this is expected and is of no consequence. 

Also, one may worry that the Wilson fermions may break parity since they have mass in 
the supercritical region. This is obviously not the case as can be seen from Figure 4 where 
the corresponding operator iJ 4 (m ) = ^p w (—mo) has no zero eigenvalues at mo = 1.9. 
Nevertheless, I measure the expectation value of \J/^y 5 \I/ for the three lattice spacings used 
in this work. The results are shown in Figure 7. Its value is of the order 10~ 6 with an 
error that is an order of magnitude larger. It is consistent with zero well within the error. 
Clearly parity is not broken and one is away from the Aoki phase |2*B] . 

The residual mass m res is measured using the ratio method (see for example ^ZT\ )• In 
Figure 8 the quenched DWF and GDWF residual masses m res vs. L s at a -1 pa 1.0 GeV 
are shown. The top points (squares) are from the quenched DWF simulations at (3 — 5.7, 
while the bottom points (diamonds) are from the quenched GDWF simulations at (3 = 4.4. 
In both cases m = 1.9 and m res was measured at mj = 0.02 (the value of m res is fairly 
insensitive to the value of m/ as can clearly be seen in Figure 11). The expected faster 
exponential decay and much smaller values of m res are evident. Already at L s = 24, m res is 
smaller by almost an order of magnitude. Furthermore, its value m res ~ 0.002 is rather nice 
for such a large lattice spacing. At L s = 40 m res pa 0.0005. The difference becomes more 
dramatic as the lattice spacing is decreased. The results for a -1 pa 1.4 GeV are shown in 
Figure 9 and for a -1 ~ 2.0 GeV are shown in Figure 10. Figure 9 is most interesting since 
it corresponds to the lattice spacing of dynamical N t = 8 thermodynamics at the critical 
temperature. At L s = 24 m res pa 0.0006 and at L s = 32 m res pa 0.0002. 

As mentioned above m res is rather insensitive torn/. This is shown in Figure 11 for the 
GDWF simulations at all three values of (3. 
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Regarding chiral symmetry, one of the most telling observables is the pion mass. At the 
end of the day any improvement method should result to small pion masses. In Figure 12 
the GDWF pion mass squared vs. rrif from quenched GDWF simulations is shown. The 
squares are the measured data points, the straight line is a least x 2 fit an d the star is the 
rrif = extrapolated point. Here f3 = 4.4 corresponding to or 1 1.0 GeV , m = 1.9 
and L s = 16. The straight line fit intersects the x-axis atm/ ~ —0.004 which is consistent 
with the value of m res ~ 0.006 from Figure 8. The data for (3 = 4.6 corresponding to 
a -1 ~ 1.4 GeV , are shown in Figure 13 and for (3 = 4.8 corresponding to a" 1 « 2.0 GeV 
are shown in Figure 14. In both cases the straight line fit gives « at mj = within 
the error bar. 

5 Light pions in rugged terrain 

In this section I present a simulation of quenched GDWF at a -1 = 1.356(75) GeV , L s = 24 
and rrif = 0.005. The space-time volume is 16 3 x 32 and m = 1.9. Think of this as the 
icing on the cake or the cherry on top of the sundae. The point is that this simulation 
should closely resemble a dynamical GDWF simulation near the N t — 8 critical point of 
QCD thermodynamics. Also notice that in order to mimic the way dynamical simulations 
are analyzed I calculated the lattice spacing using the p mass at the bare quark mass 
rrif — 0.005 and not at the rrif — extrapolated value. 

I find that m res = 0.00064(4) which is about 10% of the explicit quark mass rrif — 0.005. 
And finally the pion mass is under control. I find = 140(40) MeV at rrif = 0.005. 

6 Topological issues 

The resolution of the flow diagrams of Figure 4 does not allow me to measure the net 
index but does allow me to observe that the number and location of crossings change from 
configuration to configuration. This indicates the desired instanton anti-instanton activity. 
Measuring the net index requires more computational resources and is outside the scope of 
this paper. However, a few remarks are in order. 

If the update algorithm "smoothly" transforms the gauge field configuration then the 
net index will have to change in a smooth way too. In that case the heavy Wilson fermion 
determinant will prohibit any flow line from crossing through m and as a result the net 
index will not be able to change. The simulation will generate configurations with the 
same net index as the initial configuration and will not be able to tunnel between sectors. 
For example, for an ordered start the net topology will always remain zero. This does 
not change the ability of the simulation to generate crossings (for example as in Figure 
4). It simply means that the appearance / disappearance of an instanton will always be 
accompanied by that of an anti-instanton of some size at some location. 

Because I could not measure the net index it is not clear if the HMC Phi algorithm 
that I used in this work is capable of topologically non-smooth gauge field evolution that 
would generate tunneling between sectors. However, this is an algorithmic issue and is not 
particular to GDWF. For example, net index change is suppressed as the lattice spacing 
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gets smaller irrespectively of using or not using GDWF. The gauge action barriers between 
topological sectors are a property of QCD. This has not been identified as a problem yet 
because the couplings used in today's simulations are not weak enough. Even so, algorithms 
that are able to tunnel between sectors have already been proposed (for example see ^S] 
and references therein). 

Therefore if net index change is important to the physics at hand one should check the 
net index properties of HMC Phi under GDWF. If they are not satisfactory then one should 
use a properly augmented evolution algorithm. In any event although the subject of net 
topology is of interest in many cases one only needs to stay within sector zero provided that 
the volume is large for the physics at hand. For a very nice investigation on the subject 
see |T5j . 

7 Algorithmic and computational costs 

The algorithmic implementation of a new method is of course a matter of effort and not a 
matter of concept. However, this effort is not always small. A method that is simple and is 
easily implemented as an extension of existing methods is highly desirable. This is the case 
for GDWF. Any QCD code that employs standard DWF or related fermions already has 
an implementation of the Wilson p w (x, x'; //) operator and evolution force term. Typically 
one should simply have to add to the DWF force term the Wilson force term and to the 
corresponding energy function the extra Wilson fermion energy. 

The additional computational cost of adding 2-flavors of heavy Wilson fermions to a 
2-flavor dynamical DWF simulation is obviously at most an 1/L 8 fraction. For example, 
at L s = 24 it should be less than a ~ 5% addition to the computation cost. 

8 Curious thoughts 

GDWF bring up a few thoughts that are worthwhile pondering about. 

Clearly the addition of heavy fermions to the theory is an alteration in the ultraviolet 
regime. But in the GDWF case it dramatically affects the infrared properties of the theory. 
After all this was its purpose. This is an example of physics above the cutoff affecting 
the physics far below the cutoff. These "spectator" fermions may be pointing to some 
interesting theoretical directions. 

GDWF were described as the addition of two extra Wilson fermions with supercritical 
mass. This leads to a fermionic determinant that directly suppresses the unwanted gauge 
configurations. In this light GDWF are simply a device for better simulations of QCD on 
the lattice. Perhaps there is nothing more to it than that. But it is interesting that GDWF 
can be interpreted in two additional and different ways: 

Another way to think about GDWF is to consider the log of the two flavor Wilson 
determinant (with mass —mo) as part of the pure gauge action. As mentioned above it is 
irrelevant and therefore it is a valid addition to the lattice gauge action. Thinking about it 
this way opens the method to be applied to the existing wealth of improvement techniques 
for DWF and related fermions. 
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Yet another and perhaps most curious way to think of GDWF is to extend the five 
dimensional DWF Dirac operator to include two additional diagonal terms along the fifth 
direction. These terms can be inserted "behind" the walls. The determinant of this new 
Dirac operator is the product of the original DWF determinant times the two flavor Wilson 
determinant. More specifically, one will have L s + 2 fermion fields (but still L s Paulli- 
Villars fields) and the DWF Dirac operator of eq. ^ will have to be augmented to the 
GDWF operator by extending it along the fifth dimension to include two p terms (see eq. 
I2J) along the diagonal. In matrix form on the fifth dimension index the GDWF operator is 
(for simplicity I have set rrif = 0, all blank entries are zero and Pr } l is given in eq. 0]): 



IP 



\ 



p-l 

Pl 



Pr 



P-l 

Pr, 



Pr 
P-l 



\ 



(8) 



PJ 



In some sense this is a natural extension of DWF. One can think of this as having 
fermions beyond the walls that do not communicate directly with the fermions inside the 
walls. Their presence is felt only through their coupling to the gauge field in the bulk. 
Notice that the fifth direction mid-point reflection property of DWF is maintained this way. 
However, also notice that this way one adds two heavy Wilson fermions per DWF flavor. 
For example, 2-flavor QCD will have an additional four flavors of heavy Wilson fermions. 
In the numerical simulations of this work I only used 2-flavors of Wilson fermions. However, 
using more flavors should not present problems. 



9 Open issues 

As mentioned earlier one does not expect dynamical simulations of GDWF to change the 
quenched results by much, simply because the lattice spacings considered here are large and 
are the ones of current interest to dynamical simulations. Nevertheless, dynamical tests of 
GDWF would be very welcomed. 

The 2-flavors of Wilson fermions that were added in the action are designed to produce 
a determinant det[H^\ 2 . This suppresses the configurations for which the transfer matrix 
of DWF has eigenvalue 1. However, if 4 is the transfer matrix Hamiltonian for DWF with 
a continuous fifth direction. Here I used DWF with a discrete fifth direction. To be sure, 
the continuous and discrete transfer matrix Hamiltonians have the same zeroes. But the 
continuous transfer matrix Hamiltonian used here may be more effective for the overlap- 
Neuberger fermions ^21 EH] as well as their variants and improvements (see 121 13 lU 13 El 
EUHlin] and references therein) since it directly relates to them. On the other hand, for DWF 
with discrete fifth dimension one may want to use instead of a standard Wilson fermion 
an augmented Wilson fermion that will perhaps further improve the chiral properties of 
DWF. The exact form of H± for the discrete case is known exactlv [T2"l I17j. Both of these 



11 



observations deserve further scrutiny. Furthermore, investigating the applicability of GF 
to other related fermion regulators [2H] is of great interest. 

Throughout this work the reader may have been wondering why one should only consider 
two flavors of Wilson fermions. Why not more? This is indeed a very interesting question 
deserving further exploration for the obvious extra benefits it may provide in the cost 
of numerical simulations, but also for purely theoretical reasons. Is a theory with many 
spectator fermions above the cutoff of any interest? 

As mentioned earlier, in production simulations GDWF should be used with improved 
gauge actions (for example Iwasaki). It is not clear if this will further improve the chiral 
properties. This certainly deserves more investigation. However, it should be used anyway 
because it more closely resembles the continuum gauge action. 

As discussed in section |B]the subject of net topology change is beyond the scope of this 
work. However it is an interesting open issue that needs further investigation. Besides the 
algorithmic approach discussed in section |H1 one may want to consider alternative ideas. 
For example, it has been proposed that using lp w 2 + h 2 , where h is a small real number, 
instead of p w 2 for the two heavy Wilson flavors will produce more frequent net-topology 
changes as h is increased from zero Clearly h interpolates between GDWF (h = 0) 

and DWF (h » 1). 

Finally, simulations of the J\f = 1 Super Yang-Mills theory are possible [201 Ej arid 
have been performed [22] using DWF. Issues similar to the chral symmetry of QCD are 
present and therefore one would expect that GDWF should work there as well. Further 
exploration is needed. 

10 Conclusions 

In this work I demonstrated that the chiral properties of Domain Wall Fermions (DWF) 
in the large to intermediate lattice spacing a regime of QCD, 1 GeV ^ a -1 < 2 GeV , 
are significantly improved by adding to the action two standard Wilson fermions with 
supercritical mass equal to the negative DWF five dimensional mass m . Here I used m = 
1.9. I performed quenched DWF simulations and showed that the eigenvalue spectrum of 
the transfer matrix Hamiltonian develops a substantial gap and that the residual mass 
decreases appreciatively (by about an order of magnitude). I used this key property to 
name this fermion regulator as Gap Domain Wall Fermions (GDWF). For other DWF 
related fermions (for example overlap-Neuberger fermions) I argued that this method of 
adding "gap fermions" (GF) should also improve their chiral properties. 

In order to make sure there are no adverse "side-effects" , I confirmed that local topology 
changing (measured using the overlap index method) remains active; instantons and anti- 
instantons with size larger than a lattice spacing are active in the gauge field configurations. 
Also, I confirmed that the hadron spectrum of the added Wilson fermions is above the 
lattice cutoff and therefore is irrelevant. Furthermore, I checked that parity is not broken 
and therefore the added Wilson fermions are well outside the Aoki phase. All numerical 
simulations in this work were performed on 512- node IBM BlueGene/L supercomputers. 

Furthermore, I argued that the results of this work should also hold for dynamical 
GDWF since I considered rather large lattice spacings which are typical to currently possible 
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dynamical DWF simulations. In particular, the middle of the range of the lattice spacings 
1 GeV < a -1 < 2 GeV approximately corresponds to the N t = 8 critical coupling of the 
QCD thermal phase transition. This makes the method very appealing for dynamical QCD 
thermodynamics. 
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Appendix 

The evolution in all simulations was done using the standard HMC Phi algorithm. In all 
cases the trajectory length was set to 0.5. The step size was set to 0.005 or 0.01 depending 
on the acceptance rate. The acceptance rate in all cases was above 85%. An initial 200 
trajectories were used for thermalization. A total of 50 to 100 measurements were done 
for all propagators. All measurements were separated by 20 trajectories. The conjugate 
gradient residual was set to 10~ 7 . 

The 2-flavor Wilson simulations are obviously full dynamical simulations in the Wilson 
fermions (they are quenched only with respect to the DWF fermions). The conjugate 
gradient residual was set to 10~ 6 . The conjugate gradient iterations for the evolution 
varied between 100 and 200. 

All fits in this work including the propagator fits to extract the hadron masses (not 
shown) had a x 2 per degree of freedom less than 1. 



Tables and Figures 

In this section the numerical simulation results are presented in the following tables and 
corresponding figures. 
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p 


m f 


m p 


m N 


4.4 


0.0 


0.720(12) 


0.969(23) 


4.4 


0.02 


0.759(10) 


1.058(26) 


4.4 


0.04 


0.789(10) 


1.115(14) 


4.4 


0.06 


0.815(11) 


1.203(14) 


4.4 


0.08 


0.868(10) 


1.275(12) 


4.6 


0.0 


0.541(9) 


0.712(25) 


4.6 


0.02 


0.596(7) 


0.823(22) 


4.6 


0.04 


0.636(19) 


0.935(20) 


4.6 


0.06 


0.690(6) 


1.036(14) 


4.6 


0.08 


0.753(7) 


1.154(15) 


4.8 


0.0 


0.399(16) 


0.558(30) 


4.8 


0.02 


0.441(27) 


0.640(47) 


4.8 


0.04 


0.522(9) 


0.776(19) 


4.8 


0.06 


0.589(7) 


0.890(9) 


4.8 


0.08 


0.644(6) 


0.993(11) 



Table 1: The p (m p ) and nucleoli (mjv) masses in lattice units for quenched GDWF with 
2-flavor Wilson fermions, V = 16 3 x 32, L s = 16 and mo = 1.9. The m/ = 0.0 data are 
from linear extrapolation. Plotted in Figure 1. 



(3 DWF 


(3 GDWF 


m p DWF 


m p GDWF 


m N DWF 


m N GDWF 


a' 1 GeV 


5.7 


4.4 


0.756(22) 


0.720(12) 


1.03(6) 


0.969(23) 


1.070(18) 


5.85 


4.6 


0.549(14) 


0.541(9) 


0.74(2) 


0.712(25) 


1.423(25) 


6.0 


4.8 


0.404(8) 


0.399(16) 


0.566(21) 


0.557(30) 


1.930(77) 



Table 2: Matching the scale of quenched GDWF with the scale of quenched DWF. The m p 
and itln data are from m/ = extrapolations and are in lattice units. The a -1 is derived 
from the GDWF m p mass. The DWF data are from reference |2Z1- The parameters are as 
in tabled Plotted in Figures 2 and 3. 



a' 1 GeV 


rrin GeV 


m p GeV 


rriN GeV 


< ^75* > 


1.070(18) 


1.373(25) 


1.530(26) 


2.97(24) 


-0.59(13.77) 10~ 6 


1.423(25) 


1.791(31) 


1.906(33) 


3.80(15) 


-1.65(10.65) 10-° 


1.930(77) 


2.318(93) 


2.406(97) 


4.71(23) 


-3.16(11.94) 10" 6 



Table 3: The m w , m p and in GeV measured using 2-flavor dynamical Wilson fermions 
(both for the propagator and the sea quarks) with mass \l — —1.9 for three values of the 
lattice spacing. The a -1 is derived from the GDWF m p mass and is the same as in table 
121 All hadron masses are above the cutoff. Also in the table the expectation value of the 
Wilson fermion ^75^ condensate is shown in lattice units. It is consistent with zero well 
within the error. Plotted in Figures 6 and 7. 



16 



a- 1 = 1.070(18) GeV 


L s 


m res DWF 


m res GDWF 


8 


0.034860(885) 


0.024820(705) 


16 


0.016765(610) 


0.006180(285) 


24 


0.010395(545) 


0.002215(120) 


32 


0.007595(585) 


0.001075(45) 


40 


0.006565(270) 


0.000550(30) 



Table 4: The residual mass m res for quenched DWF and GDWF for various L s values. 
The parameters are V = 16 3 x 32, /9 DWF = 5.7, (3qdwf = 4.4, m = 1.9, and rrif = 0.02. 
Plotted in Figure 8. 



a- 1 = 1.423(25) GeV 


L s 


m res DWF 


m res GDWF 


8 


0.013650(325) 


0.012700(450) 


12 


0.006950(385) 


0.004785(115) 


16 


0.004585(155) 


0.002095(65) 


20 


0.003255(110) 


0.001035(75) 


24 


0.002660(100) 


0.000595(34) 


32 


0.001950(75) 


0.000231(14) 


40 


0.001560(65) 


0.000102(7) 


48 


0.001380(60) 


0.000051(4) 



Table 5: The residual mass m res for quenched DWF and GDWF for various L s values. 
The parameters are V = 16 3 x 32, Pbwf = 5.85, /Sqdwf = 4.6, m = 1.9, and rrif = 0.02. 
Plotted in Figure 9. 



q-1 = 1.930(77) GeV 


L s 


m res DWF 


m res GDWF 


8 


0.005810(345) 


0.006640(405) 


16 


0.001085(70) 


0.000650(30) 


24 


0.000495(35) 


0.0000960(85) 


32 


0.000310(25) 


0.00002055(195) 


40 


0.000225(25) 


0.00000525(45) 



Table 6: The residual mass m res for quenched DWF and GDWF for various L s values. 
The parameters are V = 16 3 x 32, /9 DWF = 6.0, /Sgdwf = 4.8, m = 1.9, and rrif = 0.02. 
Plotted in Figure 10. 
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771/ 


m res , (3 = 4.4 


m res , (3 = 4.6 


m res , (3 = 4.8 


0.02 


0.006180(285) 


0.002095(65) 


0.000650(30) 


0.04 


0.005665(145) 


0.001945(105) 


0.000655(25) 


0.06 


0.005720(160) 


0.001900(80) 


0.000650(15) 


0.08 


0.005665(105) 


0.001995(80) 


0.000655(15) 



Table 7: The quenched GDWF residual mass m res for various values of to/ and for three 
values of (3. The parameters are V = 16 3 x 32, L s = 16 and m = 1.9. Plotted in Figure 



TO/ 


ml, (3 = 4.4 


ml, (3 = 4.6 


ml, (3 = 4.8 


0.00 


0.021(13) 


-0.005(12) 


-0.0001(94) 


0.02 


0.137(12) 


0.086(8) 


0.063(6) 


0.04 


0.227(9) 


0.166(15) 


0.117(9) 


0.06 


0.343(12) 


0.259(14) 


0.182(13) 


0.08 


! 0.451(8) 


0.363(22) 


0.255(18) 



Table 8: The quenched GDWF pion mass squared ml for various values of to./ and for 
three values of (3. The to/ = data are from linear extrapolation. The parameters are 
V = 16 3 x 32, L s = 16 and m = 1.9. Plotted in Figures 12, 13 and 14. 
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1.5 





Figure 1. The quenched GDWF nucleoli (top points) and the rho (bottom points) mass 
in lattice units vs. rrif at L s = 16 for three values of (3. The extrapolated rrif = rho mass 
is used to set the scale. 
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4.8 



0.0 
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5.9 
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Figure 2. Matching the scale of quenched GDWF with the scale of quenched DWF. The 
nucleon (top points) and the rho (bottom points) mass in lattice units for three values of /3 
is shown. The diamonds are for DWF from reference [27] and the crosses are for GDWF. 
Also the j3 values for DWF are the bottom numbers of the x-axis while the GDWF f3 values 
are the top numbers of the x-axis. 
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Figure 3. The inverse lattice spacing in GeV vs. the quenched GDWF (3 values used n 
this work. 
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Figure 4. The ten smallest magnitude eigenvalues of H^mo) vs. m . An aggregate of the 
results from 20 independent configurations is plotted in each plot. The left column is from 
0-flavor Wilson (quenched) simulations while the right column is from 2-flavor Wilson with 
mass —itlq. Horizontally the 0-flavor and 2-flavor f3 values correspond to the same lattice 
spacing (Figure 2). From top to bottom a -1 pa 1.0 GeV , 1.4 GeV and 2.0 GeV . Notice 
the difference in the y-axis scale for the different lattice spacings. 
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Figure 5. The 100 smallest magnitude eigenvalues of #4 (mo) at mo = 1.9 were measured 
for 110 independent configurations for (3 = 5.85 and 0- Wilson flavors and for (3 = 4.6 and 
2- Wilson flavors. Both have a -1 ~ 1.4 GeV . The fraction of eigenvalues with values 
between A and A + dX is plotted vs. A. Here dX ~ 0.007. 
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1.0 1.2 1.4 1.6 1.8 2.0 

1/a (GeV) 



Figure 6. The pion (diamonds), rho (squares) and nucleon (stars) masses vs. a -1 in GeV 
measured using 2-flavor dynamical Wilson fermions (both for the propagator and the sea 
quarks) with mass \x = —1.9. The scale is set using the GDWF rho mass. The straight line 
marks the cutoff. Masses above that line are above the lattice cutoff. 
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Figure 7. The expectation value of the ^75 \l/ condensate in lattice units vs. a -1 in GeV 
using 2-flavor dynamical Wilson fermions (both for the propagator and the sea quarks) 
with mass n = —1.9. It is consistent with zero well within the error. The lattice spacing is 
set using the GDWF p mass. The straight line marks the y-axis zero line. 
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Figure 8. The quenched DWF and GDWF residual masses m res vs. L s at a~ l « 1.0 GeV . 
The top points (squares) are from quenched DWF simulations at f3 = 5.7 while the bottom 
points (diamonds) are from quenched GDWF simulations at /3 = 4.4. In both cases mo = 
1.9 and m res was measured at mj = 0.02. 
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Figure 9. The quenched DWF and GDWF residual masses m res vs. L s at a~ l m 1.4 GeV . 
The top points (squares) are from quenched DWF simulations at (5 = 5.85 while the 
bottom points (diamonds) are from quenched GDWF simulations at (3 = 4.6. In both 
cases m = 1.9 and m res was measured at m; = 0.02. 
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Figure 10. The quenched DWF and GDWF residual masses m res vs. L s at a~ l « 2.0 GeV . 
The top points (squares) are from quenched DWF simulations at j3 = 6.0 while the bottom 
points (diamonds) are from quenched GDWF simulations at (3 = 4.8. In both cases mo = 
1.9 and m res was measured at rrif = 0.02. 
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Figure 11. The quenched GDWF residual mass m res vs. m/ at L s = 16 for three values 
of p. 
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Figure 12. The pion mass squared vs. m/ from quenched GDWF simulations. The 
squares are the measured data points, the straight line is a least \ 2 fit and the star is the 
TTif = extrapolated point. Here (3 = 4.4 corresponding to a -1 ps 1.0 GeV , m = 1.9 and 
L s = 16. 
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Figure 13. The pion mass squared vs. m/ from quenched GDWF simulations. The 
squares are the measured data points, the straight line is a least \ 2 fit an d the star is the 
uif = extrapolated point. Here (3 = 4.6 corresponding to a -1 m 1.4 GeV , itlq = 1.9 and 
L s = 16. 
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Figure 14. The pion mass squared vs. m/ from quenched GDWF simulations. The 
squares are the measured data points, the straight line is a least \ 2 fit an d the star is the 
TTif = extrapolated point. Here (3 = 4.8 corresponding to a -1 ps 2.0 GeV , m = 1.9 and 
L s = 16. 
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